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We calculate the one-loop effective action of the 5(7(2) Yang-Mills theory for center- vortex con- 
figurations, both in 3d and 4ci. We find that in both cases there are minima of the effective action, 
corresponding to vortices of the transverse size approximately 4/(73 smd 1.7/A-jf7g-, respectively. The 
values of the effective actions at the minima are negative, suggesting that the Euclidian vacuum 
may be unstable with respect to creation of vortices. 
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^ ' I. INTRODUCTION 

. Recently there have been certain indications from lattice simulations that center-of-group vortices P, may be 
responsible for the area behavior of large Wilson loops, i.e. for confinement, see e. g. UflSl ^'^'^ references therein. 
■<!;;j- If long closed vortices populate the Euclidean vacuum of QCD, resembling curved and entangled spaghetti in Italian 
pasta, it might serve as a microscopic mechanism of confinement. 

(N , 

^ . From the continuum side, however, vortices invoke more questions than there are answers today. First of all, 
there are no finite-size vortex-type solutions of the classical Yang-Mills equations in a non-compact space, as the 
classical equations are dimensionless whereas the functional to be minimized is the transverse energy with the 
• dimension 1/cm^. Therefore, if one finds a solution, its dilatation to larger size will have less energy. This shows 
\ that vortex- type solutions of classical equations do not exist 6]. In a non-compact space the length scale and hence 
Cn| • the vortex radius can be only set by quantum fluctuations. Therefore, vortices can at best arise as saddle points 
I of the effective action, with quantum fluctuations about a trial vortex configuration taken into account. Second, if 
the effective action contains vortex configurations as saddle points, their energy per unit length in Sd or per unit 
surface in 4d should be either small positive or, even better, slightly negative. If it is large and positive there is no 
' , reason for the vortices to appear in the physical vacuum, in the first place. Third, the typical lengths where vortices 
bend to an angle of the order of unity must be on the average larger than the characteristic transverse size of the 
, vortex. If it is not so, the fluctuation can be hardly called a vortex. Fourth, different vortices in the vacuum or 
■ different segments of the same vortex need to be on the average repulsive so that they do not glue up. If vortices are 
. !^ ' permanently merging and splitting it is difficult to consider them as adequate degrees of freedom. A pilot study of 
the 'spaghetti vacuum' in 4c? has been performed many years ago by the Copenhagen group 0, who considered 
^ ' a 2c? lattice of parallel flux tubes and then presented arguments that the lattice should "melt" into a disordered 
phase. Ultimately, one has to build a consistent statistical mechanics of the vortex vacuum in order to see that the 
'entangles spaghetti' have a preferred vacuum free energy as compared to, say, the instanton liquid. The statisti- 
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cal mechanics of instantons is quite a difficult problem by itself 0, , but vortices are far more difficult to deal with 



In this paper, we address the first and in fact the easiest questions from the above list. For vortices to be physical 
objects and not artifacts of a regularization and gauge fixing, their typical transverse sizes must be finite in physical 
units, i.e. to be of the order of l/gf in 3c? and of the order of 1/Aqcd in 4c?. Correspondingly, the energy of a vortex 
per unit length should be of the order of g| in Zd, and the energy per unit surface should be of the order of Aqqj-, in 
4c?. The appearance of Aqcd implies the transmutation of dimensions through renormalization; it can be achieved 
only with quantum fluctuations about the vortex taken into due account. We study here the effective action with 
quantum fluctuations about a trial vortex configuration being integrated out, at the one-loop level. Although the 
one-loop approximation is not the whole story, it is the first attempt to find out if there are any hints from the 
continuum Yang-Mills theory that isolated 'thick' vortices have the right to exist. We seem to get a positive answer to 
this question. It seems that quantized-flux vortices are dynamically preferred compared with the perturbative vacuum. 



In what follows, we consider an idealized isolated infinitely long and straight vortex perpendicular to the {xy) 
plane. We choose an Ansatz for the vortex profile in the transverse plane Ap_{x,y) and expand the YM action 
into the classical part, j df^xF^j^^A) / Ag\, and the part quadratic in the fluctuations about the background A^. 
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Integration over (and over ghost fields stemming from background gauge fixing) is non-trivial, as the mere 
counting of levels as compared to the free case is not obvious for the quantized-flux vortices; the difficulty is related 
to the Aharonov-Bohm effect. We perform this integration by extending the Jost scattering theory (which 'counts' 
the levels properly) to the two-dimensional case: it results in the quantum part of the effective action being finite 
after renormalization. The full effective action is the sum of the classical and quantum parts. In fact, we study 
several functional forms of the profile and find that in each case there exists a minimum of the action with respect 
to the transverse spread of the vortex. It turns out to be about (4 — h)/g1 in 3d and about (1.6 — 1.8)/Ajgg- in 4(i. 
The values of the effective action at the minima are negative in all cases considered. It indicates that the Yang-Mills 
vacuum might be unstable in regard to the spontaneous production of vortices. 

This work has certain overlap with the earlier study by the Copenhagen group 0, who started from the 
observation by Savvidy that already a constant chromomagnetic field lowers the Yang-Mills vacuum energy. 
The Savvidy vacuum has, however, a negative fluctuation mode leading to instability. The Copenhagen group has 
included the amplitude of the negative mode (interpreted as a Higgs field but actually part of the gauge field) into the 
minimization procedure and indicated that the 2d hexagonal lattice of tubes with quantized flux leads to a further 
lowering of vacuum energy. However neither the fluctuations in the unstable mode nor in the higher modes were 
included. In contrast, we take into full account all quantum fluctuations about a single trial vortex. The full one-loop 
quantum calculation performed here shows that already isolated center vortices in SU(2) and also their embedding in 
SU(3) are energetically preferred compared with the perturbative vacuum. This result comes both in Ad and 3d (the 
latter case has not been considered in Refs. |3,[3) giving thus certain support to the disordered- vortices scenario of 
confinement in both cases. 



II. THE VORTEX ANSATZ 



By definition of a Z{N) vortex, the Wilson loop in the fundamental representation of the SU{N) gauge group, 
winding around the vortex in the transverse plane assumes values being nontrivial elements of the group center, as 
the radius of the Wilson loop tends to oo: 

Pexpi * A^dx*^ ~> ... G -^w, (1) 

J \ j 

where k — 1, . . . , N — 1. In this paper, we consider the SU{2) gauge group (although we also make calculations in 
SU{3) for additional checks). In SU{2) the only nontrivial element of the center is the minus unity 2x2 matrix. 
Taking for example a circular Wilson loop of radius p in the {xy) plane around a vortex centered a.t p — \/ + = 
we see that only the azimuthal component of the Yang-Mills field is involved in eq. (Q), — eapnaA^, where Ua 
is a unit vector in the (xy) plane. One can always choose a gauge where is independent of the azimuthal angle 
(j). Generally speaking, it implies that the radial component Ap is nonzero, however, we shall neglect this component 
as the condition is compatible with Ap = 0. If Ap is nonzero it can be reconstructed from gauge invariance by 
replacing dp dpS"-'' + f°-'^^Ap. A circular Wilson loop in the J = ^ representation lying in the transverse plane and 
surrounding the vortex center is |l3j: 

Wi{p) - ^Tr Pcxpz^pA^rdf^ - cos[7r/i(p)], ^^(p) = p^ A-{p)A-ip), (2) 
For an arbitrary representation of SU{2), labeled by spin J, one has: 

W (n) - 1 sin[(2J-|-l)7rMp)] 
^'^'^-2J+1 sinMp)] • 

If p{p) tends to an odd integer at large p, the Wilson loop Wj{p) (—1)^"'. For half-integer representations it is 
minus unity, hence /i(oo) = 2n -|- 1 is the condition that one deals with a Z{2) vortex. We shall only consider the case 
when /i(oo) = 1. For simplicity we assume that only one color component of A^ is nonzero: 

Al{p) = S^'^, /i(0) = 0, mM = 1, (4) 
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where fi{p) will be called the profile of the vortex. Eq. with all the rest components of the YM field set to zero, 
is the Ansatz for the center vortex field we are going to investigate. 



III. THE CLASSICAL ACTION OF THE VORTEX 



We take an idealized straight-line (in 3d) or straight-surface (in 4d) vortex whose dimension in the longitudinal 
direction is (d = 3, 4), while in transverse plane it has a trial profile fi{p). The classical action of the vortex is 

(see Appendix A): 



Sci^s = ^ 1 d'x{F;,f = d-^-^xii Jd'x^ {B^f 



-{d-2) 



-2 / dpp 



n 2 



= h'-f ^ / dpp 
" 92 Jo 



1 a 

pop 



class 



(5) 



We remind the reader that in id the coupling constant g| has the dimension of mass, and the theory is convergent. 
In Ad the (bare) coupling g\ is dimensionless but gets an infinite renormalization from quantum fluctuations. To the 
one-loop accuracy one can parameterize: 



8^ 11 1 

^ TV In - 

54 3 A, 



M 



QCD 



(6) 



where M is the ultra-violet cutoff in a particular regularization scheme. The quantum part of the effective action 
is logarithmically divergent; the renormalizability of the theory ensures that the dependence on the UV cutoff M is 
canceled in the full effective action. 



IV. DIMENSIONAL ANALYSIS 

Let us make a quick estimate of the 'best' transverse size of the vortex po from a simple dimensional analysis. We 
denote the action per unit longitudinal dimension of the vortex by £ such that S — L^^ '^'^ £, and call £ the transverse 
energy of the vortex. If po is the scale of the profile function p{p) [e.g. p,{p) — exp(— po//o)]j we find on dimensional 
grounds that in Ad 

A 117V M 

tciass - — -T-r^ In , (7) 

P5 247r^ Aqcd 

where A is the dimensionless functional of the trial profile. 



2 f°°dp / d4i\ 
[dp] 



A 2 rd/: 
A = 7rp„ — 

Jo P 



2 

(8) 



The renormalizability ensures that the quantum part of the effective action has exactly the same coefficient in front 
of the divergent In A/ piece, but with the opposite sign, so that InAf cancels in the sum of the classical and quantum 
parts: 

Equant = ln(A/poS), S^S = £c\ass + 'S^quant = ^7"^ In(AQCDPO-B), (9) 

Pq 247r^ Pq 247r^ 

where A^B are numerical coefficients depending on the concrete functional form of the profile p(p). £cS has the 
minimum at 

p- = ^^ ^eT = -A^ii^^, e = 2.71828... (10) 
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In 3d a similar dimensional analysis yields: 



^class — 991 Equant — 1 \^^) 

with the same coefficient A as in Ad. The sign of C is apriori unknown. However, since the theory is super- 
renormalizable one can suppose that the sign of the quantum energy is the same as in Ad at Mp —^ oo, i.e. negative 
(direct calculations below confirm the expectation) . The sum f off = ^ciass + Equant has then a minimum at 

mill _ _L fll omin _ _ 2 ^ /-, ^^ 

Notice that both in 3d and Ad we get a natural result for the transverse size of the vortex and a negative sign for its 
transverse energy at the minimum. 

In a sense we have already proven the main statements of the paper put in the Abstract and the Introduction: these 
statements follow from dimensional analysis and the rcnormalizability of the theory, provided the coefficients A, B,C 
are finite, which is not obvious beforehand. In the rest of the paper we compute these coefficients for different cases, 
and look for the best proffie function /i(p). 

V. THE QUANTUM ACTION ABOUT A VORTEX 

A. Definition of tlie quantum action 

To get the quantum action for /i(p) we integrate over quantum fluctuations about a given background field A'^{p) 
(see eq. (0}) considered to be a slowly varying field with momenta k± up to certain fc„iin- Accordingly, quantum 
fluctuations have momenta k± above fcmin and arbitrary longitudinal momenta. Writing = A^ + we expand 
F^j^{A + a) in the fluctuation field up to the second order appropriate for the 1-loop calculation. The term linear 
in vanishes due to the orthogonality of high and low momenta and also because A^ is in our case independent of 
(/> and longitudinal coordinates whereas is dependent on those coordinates. 

The quadratic form for is the standard (see e.g. jljj) 

W^t = -[D\A)rS^, - 2r'F^^M), (13) 
if one imposes the background Lorentz gauge condition, 

D;\A)al = 0, D;!'iA) = d^S^" + r'^Al. (14) 
The quantum part of the effective action is, as usually, given by the small-oscillation determinants: 

5quant[A] = ^giuoii + 5ghost = \ lndet(W^;,^) - lndet(-Z?2). (15) 
Subtraction of free determinants (with zero background field) is understood. 

The full effective action is 5cff — Sdaaa + "Sgiuon + ^ghost- Since all terms are proportional to the longitudinal 
dimensions of the vortex L^^ we shall divide all terms by this quantity; it will then become an equation for the 
appropriate contributions to the transverse energy of the vortex. 



^cff ^class ^~ ^gluoii ^gtiost ; 

where all terms are certain functionals of the trial vortex profile p{p). 



(16) 
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B. The ghost operator 

For the vortex Ansatz the ghost operator (D^)"'' takes the following simple form in the cylindrical coordinates 



1 ^ ^ ^.1 
p op op 



p- 



d(f) 



del) 



The eigenvalues of the structure constants if°-^^ in the color space are 



(17) 



(18) 



{0,1,-1} for SU(2), 

{0,0,1,-1,1/2,1/2,-1/2,-1/2} for SU(3). 

It is natural to choose a plane wave basis exp(i/c|| • a;||) for the longitudinal directions 0:3,... x^, and a polar 

basis e'xp{i'm(j))Zm\p) for the transverse directions. The eigenfunctions Zm\p) satisfy the second-order differential 
equations depending on the color polarization c— 1, . . . , {N^ — 1): 



Id 9 ^ {m — acp{p)y 
pdp^dp p^ 



.Z\${p). 



(19) 



For each color polarization c one has to solve a separate eigenvalue equation with a corresponding coefficient etc from 
eq. (|18l) . The spectrum in fcj^ is continuous with the spectral density depending on m and c. By putting the system 
into a large circular box of radius i?, one makes the spectrum discrete, and the eigenvalues fc_L can be labeled by a 
discrete number n. 

The eigenvalues of the full operator —D^ are then 



(c, m, n, /c|| ) = fc^ + /c^^ ^ 



(20) 



and one has to sum over the color polarizations labeled by c, the magnetic quantum numbers m, the radial 
quantum numbers n and integrate over the continuous spectrum in fc||. In the free case (/x — 0) the eigen- 
functions of eq. H19I) are the ordinary Bessel functions J\m\{}^ii\.TanP) (the index must be non-negative to ensure 
the regularity at the origin). In this case the eigenvalues fcoJ-mn can be determined by the zeros of the Bessel functions. 

According to eq. (|15|) , the ghost contribution to the transverse energy of the vortex is the sum of logarithms of the 
eigenvalues: 



f ghost = - ^ [ln£'^(c,TO,n,fc||) - ln^;g(c,m,n, /c||)] 



c,m, 

oc 00 

E E E 

c— 1 m— — 00 n 



(2^) 



d-2 



In fc 



)-'"( 



k\ -f A:qj_ ran 



(21) 



where k'^^^an eigenvalues of the free operator in a box, with p[p) set to zero. Color polarizations with ac = 

do not contribute since at = eq. H19|l reduces to the free one. 

We first integrate over k\\ . In the Sd case the integration is finite: 



dfc|| 



In 



' k\^kl 



\kx\-\k^i\. 



(22) 



In the 4d case integration over fcy needs to be regularized. We use the Pauli-Villars regularization scheme, so that 
henceforth Aqcd = Ap.y. = e^l^'^k-^g. We have 



d^fc|| (fc^-t-fci)(fcj| + fco\-hM2) _ 1 
(2^)2 ^ {k\^kl^){k\^k\+hP) 47r 



In 



M2 



k' fln^ 
/S(,^ I in ,,2 



(23) 
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There is a general formula for summation of any function of eigenvalues becoming continuous in the limit when 
the box radius goes to infinity, see Section 6. In the continuum limit the summation can be replaced by integration 
over the spectrum with a weight being the phase shift S{k±) of the corresponding differential equation, in this case 
eq. (113: 

V [Fikn) - F{kon)] - rdkdik)^^ (24) 

„ TT Jo dk 

(the subscript _L will be henceforth suppressed). We shall describe the method of finding the phase shifts in Section 6. 
The phase shifts d{k) depend on the 'partial wave' m and on the color polarization c. Let us introduce the accumulated 
phase shift, summing over all partial waves. 



5c{k)^2Y,ScUk)+Sco{k). (25) 

m— 1 

Since ac assumes as many "plus" values as there are "minus" ones, it is sufficient to sum over non-negative values of 
m. [Alternatively, one can take only positive values of ac but then sum over all m's, positive and negative.] Using 
the linearity of eq. 1)24(1 in S{k) one can introduce the aggregate phase shift of the ghost operator, 



Sshost{k) = Y,W- (26) 

c 

It is a common function both for 3d and Ad. However, the expression for the ghost energy is different in 3d and Ad, 
as the weights with which eigenvalues are taken are different, cf. eqs. (|22I23|) . Using the general eq. ^ we obtain 
the needed transverse energy of the vortex as due to the ghost part of the quantum action : 

1 f°° 

f ghost = - / dk6,host{k) (27) 
Jo 

4host = - / dfcfcln — (5ghost(fc). (28) 
''^ Jo 1^ 



C. The gluon operator 

We now consider the gluon operator W^^, see eq. (|13|) . Its {D'^)°^^5^^ part is essentially identical to the ghost 
operator. For the other part we find the following eigenvalues: 



c racb\ 



Eig(F;,/ 



p 



({adcolor ® {/^Alspaco), 



{0,1,-1} for 3d, 
{0,0,1,-1} for 4d, 



(29) 



where the color eigenvalues ac are the same as for the ghost operator, see eq. ((TH|l . 

Therefore, the transverse eigenvalues fcj^ of the gluon operator are found from solving the differential equation for 

the eigenfunctions Zm^\k±p) where c = 1, (iV^ — 1) labels the color polarization and A = 1, ...,d labels the space 
polarization of the gluon: 



^ ^ {■m-acp{p)f 2ac(3xp'{p) 
pdp^dp p^ p 



Zt'Hp)=klcX^mZt'Hp)- 



(30) 



According to eq. (|15|l , the gluon contribution to the transverse energy of the vortex is the sum of logarithms of the 
eigenvalues: 



'S'giuon = ^ ^ [lni;^(c, A,m,n, fc||) -ln£^Q(c, A,m,n, fc|| 



(31) 



c, A.m.n.j 
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Color polarizations with ac = cancel from this equation. In Ad there are two space polarizations with /3\ = for 
which the eigenvalue equation (|Sn|l reduces to that of the ghost operator, eq. Therefore, £ghost is completely 

canceled by that part of the gluon operator spectrum. In 3d there is only one space polarization with P\ = 0, so that 
the cancellation is not complete. 

Similar to the previous subsection, we introduce the phase shifts S{k) but now corresponding to eq. 1)30(1 : 



oo 



4iuon(^) = E E ^^^(^) = 2 E -^^A^lfc) + SMk). (32) 

c A m— 1 

Notice that, unlike the ghost case, here the aggregate phase shift depends on the number of dimensions as the number 
of gluon polarizations A is d-dependent. Using the results of the integration from the previous subsection we can 
write down the gluon contribution to the transverse energy via the aggregate phase shifts: 

■^gluon — ^TT' / '^gluonC'')' (33) 

^~ Jo 

SiLn - dkk\n-6l^l^ik). (34) 



D. Full quantum energy 



Adding up the ghost and gluon contributions to the quantum energy we get: 



Equant 



I74d 
''quant 



- I dk 

TT 



1 



W(fc) - ^4?uon(fc) 



— I dk k \n —— 



1 



W(fc)-2'5giuo„(fc) 



(35) 
(36) 



The combination of ghost and gluon phase shifts in the square brackets will be called the full phase shift; for various 
cases they are plotted in Figs. 2,4 and 5. 



E. Scaling properties 

Let as assume that we have calculated the ghost and gluon phase shifts for a certain profile function fj,(p/po) with 
Po = I in some arbitrarily chosen units. If we now change po all quantities have to scale in accordance with the 
dimensional analysis of Section 4. Equations for Equant of the previous subsection can be translated into the following 
expressions for the constants A,B,C introduced in Section 4: 



A 



C = 



A 



InS 



^pI /°"^(M'(p))^ 

h P 



^ I dk 



^4?uon(fc) - '5ghost(fc) 



24 

lliV 



dk k 



241n(fc)-W(fc) 



1 24 
AllN 



dk k ln(fc/9o) 



Con(fc)-W(fc) 



(37) 

(38) 
(39) 
(40) 



As discussed in Section 4, the renormalizability of Yang-Mills theory in 4d requires that A = A: if that is satisfied 
the effective action does not depend on the ultraviolet cutoff M. This relation is, therefore, an important check of 
the numerics. It can be said that the calculation of A is a particular way of getting numerically the '11/3' of the 
asymptotic freedom law. Further on, we verify the scaling relations (|9I11|I explicitly by changing the scale po of the 
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A 

iVc = 2: 

A/ A 

B 

C 





3: 


A(iVe 




A(iVc 








B(JVc 


=3) 


C{iVe 


=2) 



C{iVc=3) 



= cxp(-po/p - p/po) 



po = 1 



0.2424 



1.0603 
1.1686 
0.1094 



1.0000 
1.0082 
0.6774 



Po 



0.2424 



1.0066 
1.0692 
0.1089 



0.9999 
1.0075 
0.6766 



P(p) = exp(-(po/p)^) 



po = 1 



2.2333 



1.0007 
1.0388 
1.1651 



1.0293 
0.8767 
0.6484 



po 



= V2 



2.2333 



0.9862 
1.0204 
1.1513 



1.0331 
0.8750 
0.6511 



P{P) = (p/Pof/iip/Pof + 1) 



po = 1 



3.3767 



1.0155 
0.8854 
1.5025 



1.0107 
0.8751 
0.6388 



po 



= V2 



3.3767 



1.0026 
0.8923 
1.4693 



1.0193 
0.8747 
0.6387 



TABLE I: Results for flux-0 and flux-1 profiles used. The table shows two columns for each profile one for po = 1 and one for 
po = 72. 



trial functions and checking that A, B, C are po-independent (see also Fig. 

Another powerful check of our numerical performance is provided by the number of colors N: the quantities A and 
B should be independent of N while C must be proportional to N. The calculations successfully pass these tests for 
the A and C within small errors, whereas the sensitivity of B on the numerics is considerably larger, see Tabled 



VI. THE PHASE SHIFT METHOD 



The Z{2) vortex field fi{p) with fj,{0) — and ^(oo) = 1 is very specific from the point of view of differential 
eqs. (|19I30|) . Near the origin p « the equations resemble the free equation for the Bessel functions J,„(A;p) while at 
infinity p oo they resemble the equation for Jm±i{kp), with a shifted index. This is the essence of quantized Z(2) 
vortices, see in this respect ref. 

This fundamental property of Z(2) vortices causes certain difficulties in finding the continuous spectrum of 
eigenvalues in a given m sector. We present here an economical method of dealing with this problem; it enables one to 
find the spectrum unambiguously. We have borrowed the idea of the method from Jost scattering theory, generalizing 
it to problems with cylindrical symmetry. The application of the phase-shift method to 1-loop calculations of the 
effective action has been, in general terms, developed in Refs. [ifil 

In the following we consider a general differential equation 



Z{p) - k'Z{p), (41) 

the equation for ghosts (|19|l and for gluons (|3U|I being exactly of that form. 

One looks for the solution which is regular at the origin and follows the oscillatory asymptotics of the Bessel function 



1 d d m'^ ^ , , 
p op op p'^ 
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at infinity: 



n,, \ ; 2 r, TT rrni . 



(42) 



This serves a definition of the phase shift S{k). If the system is put in the box of large radius R the zeros of (|42f) 
define the spectrum: 



, , S{k^) TT / 3 TO 



(43) 



where fco„ is the spectrum with zero potential. A sum of any function of the eigenvalues kn in the continuum limit 
(i? oo) becomes: 



1 1 f°° 
y2[F{k^) - F{kon)] = -^y^F'ikn)Sikn) ^ / dkF'{k)8{k). 



(44) 



This is eq. H24(l used in the previous section. 



In the free case [U — 0) the differential eq. H41|l have regular solutions Jm{kp) (the ordinary Bessel functions), as 
well as singular solutions Ym{kp) (v. Neumann functions). From those one combines Hankel functions, 

H^{kp) = J™(fcp) ± iY^{kp). (45) 
Taking into account that J,„ and Ym are both solutions of the same differential equation, their Wronskian is a constant: 



W{Jm,Yr^) = kp[J^{kp)Y:^{kp) ~ Jl^{kp)Y^{kp)] = -. 

TT 

Using H46|l one can easily verify that the needed solution Z„i{kp) of eq. H41I) satisfies the integral equation: 



(46) 



Zm{kp) = Jrn{kp)-- dp p [J„^{k p)Y„^{k p) - J„^{kp)Ym{k p)]U {p)Zm{kp) 
^ Jo 

\-^- j' dp'p'Hl-\kp')U{p')ZM) 

p 



= Im|zi/i+)(M 
= Re|if(+)(M 



1 - - dp'p'H^-\kp')U{p')Z^{kp') 

It enables one to find Zm{kp) iteratively knowing the function at smaller values of p, starting from p = 0. 
At p ^ cxD one obtains the asymptotics: 



(47) 



A cos 



TT TOTT 

kp-^ ^ + 



Re ( g«('=P-^™/2-'r/4) 



l-yj^ dp'p'H'i-\kp')U{p')ZM) 



(48) 



Putting kp — TO7r/2 — 7r/4 = 2m: or kp — TO7r/2 — 7r/4 = (2n + l/2)7r with 7i being a large integer and comparing the 
expressions we obtain the equation for the phase shift: 



tan (5(A;) 



Im/„(fc) 
Re/,n(fc)' 



U{k) = 1 - y dppHi;Hkp)u{p)z^{kp) 



(49) 



We shall refer to eq. (|^ as Jost equation as it is closely related to the theory of Jost functions developed in scattering 
theory, see e.g. p^ . 
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Check of the Jost equation 




i , I I I , I I I , 

4 8 12 16 20 

m 



FIG. 1: Check of eq. 14911 . For k — 2, a = ±1 and /3 = 1 we show the phase shift as function of m. The sohd and dashed 
hne present the phase shifts as obtained by extrapolating the difference of the zeros, while the open and filled circles show the 
result as obtained from eq. 1491 . One finds perfect agreement between the two methods. It is furthermore seen that the phase 
shifts decay rapidly in m. 

We performed a numerical check of eq. 1)49(1 by taking, as an example, the potential arising from the gluon operator 
(see eq. JSOJ), 



U{p) 



:[[ni ~ a^{p)Y — TV? — 2a[3pj-i! {p)]^ with = exp(— p — l/p). 



(50) 



We computed the phase shifts directly from the shift of zeros of Zm{kp) with respect to those of the free Jm{kp), and 
then compared with phase shifts obtained from eq. (|49|l . The comparison of the two methods is shown in Fig. ^ the 
agreement is excellent. It is also seen that Sm{k) decays fast as function of m. 

One can simplify the method of finding the phase shifts even further. The second integration in eq. (|49|l is in fact 
unnecessary. Indeed, from eq. (|47|l one has the identity: 



Zni{kp) 



Hl+\kp)f^{k) + H^^\kp)rrn{k) 



(51) 



Putting it in the Wronskian (see the definition (|46|) ') we get: 



VF(Z,„, i/(-)) = i/„(fc) I^(ijW, i?^^)) = i/„(fc) I^(r„, J„) - --f,n{k). (52) 

/ TT 

Therefore the needed phase shift can be found directly from the Wronskian composed from the solution Zm{kp) and 
the Hankel function: 



<5„(A:) Arg(/„(fc)) = Arg 



(53) 



Here Arg(re"^) = (5 is the argument of a complex number. Notice that eq. (|53|) is independent of p. Therefore, 
evaluating the Wronskian in eq. (|52|l at different values of p is a valuable check of the stability of the numerics. In 
fact, this method proves to be quite stable; we actually use it to calculate the phase shifts. 
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accumulated phaseshifts: //(p)=exp( — l/p— p) 
0.6 




4.0 



FIG. 2: Full phase shifts 5 = i5, 



2 'Jgluon 



"ghost 



for the profile fi{p) — exp(— 1/p — p): Sj^J^ik) (filled circles), 5\^J,^{k) (filled 



triangles), ^^IjC^) (open squares), 5'}^'^^{k) (open diamonds). 



VII. RESULTS 



A. Flux-O vortices 



We first consider not a Z{2) but a zero-flux vortex with /i(0) = m(oo) — 0. For definiteness we take a smooth profile 
of the form: 



H{p) = exp 



P 



(54) 



Let us briefly describe the numerical procedure. We start by defining a range R for a given value of k and m up to 
which the differential equations have to be solved. To make sure that we are working in an asymptotic regime we 
choose -Rmax = max(60.0, x\jn\i5)/k where x\jn\n is the n-th zero of the Bessel function J|„|(x). We have checked that 
for this choice the phase shifts calculated from eq. (|49|l and from eq. 153|l do agree so that we are indeed reaching the 
asymptotic regime. 

The initial condition for finding a regular solution Zm{p) is the free Bessel function which behaves as Jm{kp) « 
(fcp/2)™/r(m + 1) at small p: it is very small at large m. Therefore, for small p (actually starting from some finite 
but small p,„in — 0.00001/fc and for m > 4) we solve, instead of eq. (|41|) . the corresponding equation for the function 
B(p) = Z(p)/((p/2r/r(m+l)), 



B"{p) = -B\p) 



2m + 1 B{p){{m - ap{p)Y - m? - 2al3pp'{p)) 



(55) 



with the initial conditions 



B{Pn 



^'(Pmin) 



kp„ 



2{m+ 1) 



(^Pmin)^ 

4(m + 2) 



(56) 



This differential equation is solved up to the point Rp = 0.9a;|m|i/fc where a;|m|i is the first zero of the corresponding 
Bessel function. Starting from this point the original differential eq. (|41|l is solved, with the initial conditions: 
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FIG. 3: Scaling behavior of the effective action in three and four dimensions; The figures show the effective action in three 
and four dimensions for the profile /i(/9) = exp(— po/p — p/ Po) as a function of po- The solid line shows the solution for Nc = 2 
and the dashed line for A^c = 3. The dimensionfull constants have been set to 1, i.e. (73 = 1 in three dimensions and Aqcd = 1 
in four dimensions. The filled circle and the square come from an independent calculation with po ~ explicitly plugged 
in. The crosses denote the place of the minimum of the effective action averaged over the results for po = 1 and po — \/2, c.f. 
Table El 
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B'{Rp) nRp\"' kB{Rp) /kRpy^ 
r(m+ 1) \~) 2r(m) \~ J 

We have checked that this procedure reproduces the usual Bessel functions in the free case, when — 0. 

Having found the regular solution Zm{p), we use it to calculate phase shifts in three different ways. The first 
method is to use eq. H49f) . where the integration goes up to i?max- This procedure becomes sometimes unstable 
because for large p the integrand is oscillating. The two other methods exploit eq. (|^ . We evaluate the Wronskian 
either at 0.8i?niax or at 0.9i?niax- We have checked that except in cases where the phase shift is consistent with zero 
these three different methods agree with each other. The phase shifts actually used in the results are the ones taken 
from the Wronskian evaluated at 0.8i?niax- 

The individual phase shifts at given k are then accumulated over m, c and A. Fig. |21 shows the total aggregate 
phase shifts for N=2,3 in three and four dimensions. It is seen that in all cases 6{k) goes to zero quite fast for large 
k. For very small k values 6{k) rises again after a maximum at A; = 0.4 . 

The results are shown in Table |l| The A coefficients of the classical and the quantum parts are equal within a 
small error due to the numerics, as it should be in order to ensure that in four dimensions the results are independent 
of the UV cutoff M. It is also seen that the A coefficients are independent of the number of colors. 

As a final step we check the scaling behavior predicted by eqs. (|9I11|I . To that end we perform the same calculation 
again but for a profile /Lt(/o) = exp(— v^/p— p/V^) and compare the result with the one obtained from the profile 
p{p) = exp(— 1//9 — p) and assuming the scaling eqs. (|9I11|I . Tab. shows indeed that in 3d and M the scaling law 
is fulfilled, which is yet another important cross check of our calculation. Furthermore, it is seen that in all cases we 
find a minimum of the effective action at negative values, which indicates that the perturbative vacuum is unstable 
against the production of zero-flux vortices. The values for the minima of the scaling factor are given in Table ITU 

B. Flux-l vortices 

For the calculations with the Z{2) unit-flux vortices we use the following trick which enables us to use the same 
technique as used for the flux-0 vortices. Instead of the actual flux-l profile p{p) with /i(0) — and /i(oo) = 1, we 
consider the profile 1 — p{p). From the energy point of view nothing is changed since the classical part depends 
only on the square of the first derivative of p{p) while in the quantum part it corresponds to an overall shift of the 
summation in m. We have checked that in case of SU(2) this 'reflection' indeed produces the same results as if one 
makes the calculation directly with the flux-l profile p,{p). For the SU{2) vortex embedding into SU(3) this is not 
the case because here we have a contribution corresponding to the eigenvalues 1/2 in color space and the shift in 
integer units in the summation over m does not work. 

The numerics for the flux-l vortices is more difficult than for the fiux-0 ones. The first problem arises due to the 
fact that the phase shift exceeds tt in magnitude. As the arcus tangent is a multivalued function one then has to 
reconstruct the true function S{k) by adding suitable multiples of 27r so that the resulting function is continuous in 
k. Furthermore, general constants in terms of suitable multiples of 7r/4 have to be added so that the aggregate phase 
shifts match each other for fc ^ oo at a value close to zero. 

The second technical problem is that, literally speaking, d{k) does not go to zero with k — * oo, but rather approaches 
a small constant, see Figs. However, it seems to be a numerical artefact since, as one lowers pmin in eq. (|^ . the 

small constant becomes even smaller. Unfortunately, the numerical precision does not allow a complete annihilation 
of this spurious behavior. In fact, we have to make a cut at the point in k when the function becomes a (tiny) 
constant. The cross checks of renormalizability and the N dependence are then fulfilled. We obtain for two different 
fiux-1 profiles. 



Z{Rp) = 
Z'{Rp) = 



p<^'\p) = exp(-l/p3), 
P^'\p) = pV{p' + l), 



(57) 
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P(p) = exp(~po/p - p/po) 


m(p) = exp(-(po/p)^) 


m(p) = (P/P0)V((P/P0)' + 1) 


Nc = 2: 








[ai] 


4.441± 0.014 
-0.012± 0.000 


3.857± 0.032 
-0.150± 0.003 


4.546± 0.072 
-0.163± 0.005 


Ud) TA-ll 


1.477 ± 0.093 
-0.005 ± 0.001 


1.601 ± 0.020 
-0.040 ± 0.001 


1.855 ±0.010 
-0.046 ± 0.000 


iVc = 3: 








POmin [53"1 
c(3d) r 2i 
'^min im 


3.007± 0.007 

n no7-i- n nnn 
-U.UZ / zt u.uuu 


2.506± 0.028 

-U.OODzt U.UUo 


2.903± 0.046 
n /ini -1- n ni Q 

-U.4Ulzt U.UiO 




1.488 ± 0.093 
-0.008 ± 0.001 


1.403 ± 0.016 
-0.079 ±0.002 


1.623 ± 0.009 
-0.089 ± 0.001 



TABLE 11: Summary of the physical results with error bars. The physical result is the mean value from the two calculations 
at po = 1, \/2, while the error is given by the difference of these two values divided by the square root of two. 



the values of the coefficients in the effective action as given in TablelHand Table|n| The cross checks of renormalizabiUty 
and N-dependence are fulfilled up to 5% as seen from TableHl A gain in all cases the effective action exhibits a minimum. 
In all cases the minima are at negative values of the effective action, which means that the perturbative YM vacuum 
is unstable against the production of center vortices in the case of SU(2) in three and four dimensions and that it 
is even unstable with respect to the production of the embedding of center vortices in SU(3). Comparing the two 
profiles used it is seen that the power-like profile /x^^^ gives a slightly deeper minima than the exponential-like profile 



VIII. SUMMARY AND OUTLOOK 



We have calculated the energy of quantum fluctuations (gluon and ghost) about vortices of various profiles, both 
in 4 and 3 Euclidean dimensions. Center vortices are peculiar in that the level counting for quantum fluctuations 
as compared to the free case is not apparent. Because of this, the mere finiteness of the quantum energy in 
the center-vortex background is not evident before an accurate calculation is performed. We have developed an 
economical method of 'level counting' through phase shifts of the appropriate differential operators, which is simi- 
lar to that used in Jost's scattering theory, and expressed the energy of quantum fluctuations through the phase shifts. 

Our results surpass many consistency tests. In Ad we actually manage to calculate numerically and with a good 
accuracy the "11/3" of the asymptotic freedom law; as a result we get the 'transmutation of dimensions' and express 
the vortex transverse size and transverse energy in terms of Ajg^, see Table n The appearance of the "11/3" is 
actually a powerful test of numerical precision. Explicit calculations for different transverse sizes confirm the natural 
behavior of quantum energy, which follows from dimensions. In addition we verify that the quantum energy is 
proportional to as it should be. In our approach, it implies non-trivial sum rules for the phase shifts for gluons 
and ghosts with different color polarizations. 

The effective transverse energy of a vortex (classical plus quantum) has a minimum lower than the perturbative 
vacuum both in Ad and 3d. In Ad this result is in a qualitative agreement with the well-known Savvidy's logarithm, 
(lliVc/967r^)_ff^ ln(i?/A^ const.), for the energy of a constant chromomagnetic field H, following from asymptotic 
freedom J^J. However, the constant chromomagnetic background has a negative mode Jl, 0, so that Savvidy's 
vacuum is unstable. [Because of this negative mode the fluctuation determinant is not well-deflned, and the constant 
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FIG. 4: Full phase shifts 5 = |(5giuon — <5ghost for the profile /x(p) = exp(— 1/r^): S^^%ik) (filled circles), 5'"^Z^^{k) (filled 
triangles), 5^^^2(^) (open squares), S^^'^^W (open diamonds). The points in grey have been excluded from the analysis as the 
function 5{k) goes over to a spurious constant behavior. 



6 / 6 

accumulated phaseshifts: ij,{p)=p /{p +1) 




iiii|i<ftitiiiii 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

k 



FIG. 5: Full phase shifts 5 — ^5giuon — <5ghost for the profile fi{p) — p^' / (p^ + l): (5jY=2(^) (fiUcd circles), 3^j^'/^._^{k) (filled triangles) 



c{3d) 



^N=2i^) (open squares), d^^'^^ik) (open diamonds). The points in grey have been excluded from the analysis as the function 
6{k) goes over to a spurious constant behavior. 



)r{4<i) 
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in the argument of the logarithm can not be determined.] In contrast, the vortex has a fast-varying magnetic field, 
the fluctuation determinant has no negative mode and can be computed. Qualitatively, the energy gain is expected 
from dimensional analysis and asymptotic freedom (see section 4) but its value is a priori unknown. In 3d one does 
not have the argument of the asymptotic freedom formula, in favor of the energy gain. Therefore, the fact that we 
find an energy gain for vortices both in 3d and Ad is non-trivial. The center (flux-1) vortices give a larger energy 
gain than the non-topological flux-0 ones and from that point of view are dynamically preferred. This was also 
the conclusion of the Copenhagen group 9] who studied semi-classically a 2d hexagonal lattice made of parallel 
flux tubes. However, from our calculation we cannot exclude that fluxes of other magnitude, e.g. flux-2 (which are 
topologically equivalent to flux-0 ones) are even more preferred. 

Despite different shapes of vortex profiles considered (one exponential, the other power-like) we get close numbers 
for the energy gain. It indicates that vortices are 'soft' with respect to their deformations in the transverse 
plane. A further indication in favor of this interpretation is the fast rise of the aggregate phase at small mo- 
menta, see Figs. I2I4I5I although zero modes in a strict sense do not seem to exist, see Appendix B. That makes 
the search for the absolute minimum in the 1-loop approximation rather problematic: it will be very flat any- 
how. Simultaneously, it makes the naive interpretation of thick vortices as neatly shaped spaghetti even more difficult. 

Assuming that the vacuum is unstable with respect to the spontaneous production of individual center vortices, 
their creation can be, in principle, stabilized by i) kinetic energy of vortex bending, ii) repulsive interaction of 
neighbor vortices and/or of different segments of the same vortex. Concerning vortex interactions, a Z{N) vortex 
in a regular gauge necessarily has a nonzero azimuthal component of the Yang-Mills field, decaying as 1/ p at large 
distances from the center in the transverse plane. Two such vortices will then have unacceptable strong repulsion 
unless they are 'gauge-combed' in such a way that Arf, is zero at large distances. In such a gauge, however, one 
necessarily gets a Dirac sheet (in id) or a Dirac 3-volume (in Ad) of gauge singularities. Dirac singularities from 
two vortices intersect along a line in 3d or along a surface in Ad. An intersection of two gauge singularities in 
non-Abelian theory is, generally speaking, not a pure gauge but a physical singularity. It is not clear how to avoid 
those infinities. In any case, vortex interactions are far from being an easy topic, but they need to be clarified 
before the statistical mechanics of vortices can tell us that the 'entangled spaghetti' vacuum has a preferred free energy. 



APPENDIX A: CONVERSION INTO CYLINDRICAL COORDINATES 

The conversion into cylindrical coordinates is given by the transformations: 

a; = pcos0 y = psa\<f> (Al) 



and 



One finds: 



Ax = cos 4>Ap — sin c 
Ay ^ sin (j)Ap + cos (jyA^ (A2) 

dx = cos (pa p ~"'t> 

dy = sin (f)dp + '^'^^ ^ dff, . (A3) 



Fi2 = d^Ay - dyAx 



cos (j}dp — ^^'^ dcj, ) (sin cpAp -f cos 4>Atj,) 



COS (p 

sm(j)dp H ) {cos(l)Ap — sin 0^0) 



P 

Aa, BAa, 



(A4) 



P dp 

using Ap = 0, and with — p,{p)/ p one obtains: 

F,2^p\p)Ip ^ {Fp^f ^2[ii'{p)/pf . (A5) 
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APPENDIX B: EXCLUDING POSSIBLE BOUND STATES AT it = 



The spectrum of the ghost and gluon operators discussed above is incomplete if there exist bound states at k± — 0. 
Were the vortices an exact solution of the classical equation of motion, zero modes would be inevitable. However, 
our vortices are minima of the effective action and not of the classical one, therefore there are no special reasons for 
zero modes of the small-oscillation operators. 



To check whether such states exist we substitute p 
differential equation: 



m — 



t = 1/p in the gluon operator and obtain the reflected 

m'(iA)" 



t 



Z{l/t) = 



(Bl) 



a and P are two constants depending on the color and space polarizations of the state under consideration. Using the 
notation B(t) = Z(l/t) we take the boundary conditions for a bound state -6(0) — 0. This tells only that B has to be 
regular at the origin. Now we can check the asymptotic of the equation using an explicit flux-1 profile, for example, 
= exp(— 1/p"^). One obtains asymptotically for large p and correspondingly for small t: 



ld_ d_ 



(m — ay 



Boit) = 



For this equation we require power series in t. The only non-trivial regular solution to this is Bo{t) 
Therefore, we use the ansatz: 

B{t) = . 
One obtains the differential equation for the amplitude: 



Id d \m-a\d {m - ap{l/t)f - [m - af + 2a(3p' {\/t) /t 



A{t) = 



(B2) 

(B3) 
(B4) 



The boundary conditions are now A{0) — 1, ^'(0) = to make sure that the asymptotical behavior is correct. Explicit 
calculation shows that these states behave as: 



p I™! for small p „ „ 
-Im-Qi r 1 tor m > 



(B5) 



The behavior at small p is clear because of asymptotic arguments like the one used for small t in eq. (|B2p . The 
explicit calculation shows that only irregular solutions exist for m > 0. In the case m = the solution is at least 
logarithmically divergent for large p, as |a| — 1, 1/2. Therefore, there are no localized solutions with zero fcj_. 
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